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' Abstract 

We study the dynamics of many charges interacting with the Maxwell field. The 
particles are modeled by means of non-negative distribution functions /"*" and /~ 
representing two species of charged matter with positive and negative charge, re- 
. spectively. If their initial velocities are small compared to the speed of light, c, 

I then in lowest order, the Newtonian or classical limit, their motion is governed by 

' the Vlasov-Poisson system. We investigate higher order corrections with an explicit 

■ control on the error terms. The Darwin order correction, order jv/cp, has been 

proved previously. In this contribution we obtain the dissipative corrections due to 
radiation damping, which are of order |z;/c|^ relative to the Newtonian limit. If all 



9" 



X 



particles have the same charge-to-mass ratio, the dissipation would vanish at that 
order. 



1 Introduction 

In classical electrodynamics it is well known that accelerated charges loose energy by radi- 
ation and there is a large amount of literature concerning effective equations which include 
effects due to radiation damping without giving a completely relativistic description of 
the system of fields and charges. 

A similar but more involved situation occurs in the theory of general relativity where 
accelerated matter emits gravitational radiation and is thus damped. The probably best 
studied example is the Hulse- Taylor binary pulsar consisting of a strongly self-gravitating 
system of two stars rotating about their common center of mass. Due to the difficulties 
as the non-linearity and the necessity of finding appropriate coodinates it seems to be out 
of reach to treat a system like the one already mentioned within the full theory. Hence, 
it is desirable to have effective equations valid in certain limits as in the electromagnetic 
case. 

In many applications, as e.g. the Hulse- Taylor binary pulsar, the occurring velocities 
are small compared to the speed of light. Thus, it is a natural strategy to expand the 
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metric in powers of |f /c|, (c denotes the velocity of light). The contribution in order zero 
corresponds to the non-relativistic limit where gravity is governed by Newtonian theory. 
Therefore, higher order corrections are usually addressed as post-Newtonian approxima- 
tions. For an overview concerning post-Newtonian expansions see \J ,. Whereas it is 
relatively straightforward to establish relations between the full system and the equations 
of approximation, it is much more difficult to give the relation between the solutions of 
the two sets of equations. While order zero is done in for asymptotically flat solutions 
any further progress seems to be difficult at this point. 

For this reason it seems to be useful to investigate the very similar but less involved 
system of charged matter coupled to electromagnetic fields. In [T7j and [TH| the first post- 
Newtonian approximations of the Abraham model, a model consisting of single charged 
particles coupled to the Maxwell fields which they create collectively, are considered yield- 
ing the Darwin corrections, order |f/cp, and radiation corrections, order |f/cp with re- 
spect to the Newtonian limit. Explicit estimates of the error terms are given. 

In the present paper we choose a model of many particles governed by a statistical 
approach. For sake of simplicity we assume that there are only two different species of 
matter with mass normalized to unity and charge normalized to plus unity and minus 
unity, respectively. These distributions of the large number of particles in phase-space 
are modeled through the non-negative distribution functions and /~, = 
depending on time t G M, on position x G and momentum p G M^. The dynamics is 
governed by the relativistic Vlasov- Maxwell system. 



dtf^ + P ■ V,/± ±{E + c-'p X B) ■ Vp/± = 0, 
cV X E = -dtB, cV X B = dtE + inj 

V-E = ATcp, V-B = 0, 

P ■■= fir-ndp, J := fpir-. 



dp 



(RVM, 



Here 



p= (l + c-V)-^/^pGM^ (1.1) 

is the relativistic velocity associated to p. The Lorentz force E + c^^p x B realizes the 
coupling of the Maxwell fields E[t,x) G and B{t,x) G to the Vlasov equation, and 
conversely the density functions enter the field equations via the scalar charge density 
p{t,x) and the current density j{t,x) G M^, which act as source terms for the Maxwell 
equations. The parameter c gives the speed of light for given units of time and space of 
the physical system represented. As usual we shall deal with the limit of small velocities 
by letting c ^ oo. Some background on this procedure is given Section |21 In order to 

one has to prescribe initial data for the densities and 
-- 0, 



1-1/2, 



give the Cauchy problem of JRVMt 
the fields at a certain time t, say t 



f^{0,x,p) = rc'^{x,p), E{0,x)=E:{x), B{0,x)=B:{x). 



1.2) 



Henceforward in our notation we will only express the initial data's dependency upon the 
light velocity by the subscript c while the dependency of the solution will be suppressed. 

In the next section we shall describe the post-Newtonian expansions used in this paper. 
First we define a "naive" post-Newtonian expansion of ( RVMe ) , see fj2.4|) - ()2.6|) : whereas 
this expansion is well defined for all orders in and the solutions of the expansion 
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equations up to order are good approximations of the solution of the Cauchy problem 
of ( RVMe ) if the initial data are well adapted, see and ^T] for the Newtonian limit 
and [5^ and ]JT^ for the corrections in order c~^, it does not include damping effects due 
to radiation which occur in the order c~^, see jH]. Hence, we introduce a more sophisti- 
cated expansion up to order containing a radiation reaction term, see ()2.13|) below. 
Assuming only one species of matter, say /~ = 0, this term would vanish reflecting that 
there is no radiation in this order of in that case. As in the case of single charges one 
has to circumvent the occurence of so-called run-away solutions in the resulting dynam- 
ics. Details are explained in Section [2.1.11 We show that the resulting effective dynamics 
is well defined, at least locally in time, and we give some nice further properties of the 
solutions. Proposition 12.11 and Lemma 12.21 Furthermore, we give an explicit expression 
of the fields as function of the values of the densities; in this expansion the fields are not 
a degree of freedom anymore but are enslaved by the densities, see formulas ()2.25a|l and 
(lOSbl) . 

In Section El we state our main results concerning the comparison of solutions of our 
expansion and solutions of JRVMc ) . 

The phase-space of the comparison dynamics defined in Section |21 consists of a pair 
/°'^, of smooth non-negative functions with compact support defined on x M^. The 
field quantities are to be computed from this densities e.g. by means of the formulas 
()2.25ap and ()2.25b|) . On the other hand treating the Cauchy problem of ( RVMg ) we also 
have to specify the initial fields E° and B° and it is the question for which choice of 
initial data the dynamics of Section |2l is a good comparison dynamics and for which not. 
Surely, it is possible to choose initial data for the fields such that the densities of the 
two dynamics evolve in completely different ways. In Section ITTl we will adapt the initial 
fields for the Maxwell dynamics from the comparison dynamics, see formula pC|) . That 
means for given data for the particle densities we compute special fields by means of the 
formulas ()2.25a|) and ()2.25bjl and impose these fields upon the initial values of the Maxwell 
fields. That has the advantage that, from a mathematical point of view, existence and 
uniqueness theorems of local in times solutions for both dynamics are at hand, see jTH!, 
[21j, 0- We prove that the error between the solutions of the two systems is of order 
0(c~^), see Theorem 13.21 We want to mention that in [21], JI] and j3] the fields are 
adapted up to the relevant orders in the same way. 

There are two drawbacks of this method. Post-Newtonian expansion is in essence 
an expansion of the relativistic velocity p and the retardet time t — c~^\x — y\- It is 
clear that assuming localized sources the expansion of the retardet time is only a good 
approximation in the near zone of the source where |a; — ?/| <C c. This is reflected in the 
fact that the estimates of the fields in Theorem 13.21 and 13.41 are only local in the space 
variable x. For this reason also the adapted initial fields are only reliable in the near 
zone, in particular they are not of finite energy. From a more physical point of view it 
also seem to be questionable to use the Cauchy problem at all. Recalling the motivation 
of post-Newtonian expansions one is more interested in localized systems, isolated from 
the rest of the world, which have already evolved for a long time with small velocities. 
Therefore the Cauchy problem might not be the right formulation since it is not clear 
how to incorporate these properties into the initial fields. In physics textbooks isolated 
systems are characterized by the absence of incoming radiation, that is energy coming 
into the system from past null infinity by means of electromagnetic fields, for a rigorous 
definition see [5]. For given sources fields free of incoming radiation are usually calculated 
by means of the retardet potential. Past null infinity is that region of space-time which 
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is reached in the direction of backward hghtcones. In Section 13.21 we consider a famihy 

passing through /°'^ at time t = where in 
the electromagnetic fields are computed by 



of solutions, parametrized by c of JRVM; 



contrast to the Cauchy problem of JRVM, 



means of the retardet potentials alone, see jretRVM, 



Because it is our goal to model 
slow systems we assume that the momenta are bounded uniformly in c > 1 and time 
t G M, see Assumption I3.3r b). It is not the aim of this paper to investigate existence of 
such solutions with mathematical rigor, instead we will just assume their existence and 
some nice properties used in the sequel, see Assumption 13 .31 Note however that in [Oj the 
existence of global solutions is proved for small /°'^. Furhermore it is shown that such a 
solution is also unique in a certain class of "nearly free streaming solutions" and is free of 
incoming radiation. We also want to emphasize that the Larmor formula has been proved 
for this system, see (12.1111 and [0] . The underlying physical picture is that in the absence 
of incoming radiation from outside the system any solution of ( RVMe I will approach a 
solution of ( ret RVMe ). That means that solutions of (rctKVMc' constitute a kind of initial 
layer. 

We prove that the error between solutions of our comparsion dynamics and solutions 
of dretRVMcT ) is of ordcr c~^, see Theorem 13.41 

In Section 13.41 we collect some more notations used in the proofs of Theorem 13.21 and 
Theorem 13.41 The proof of Theorem 13.21 is elaborated in Section 0] while the somehow 
cumbersome computations of some representation formulas for the fields is outsourced to 
the Appendix IHl The proof of Theorem 13.41 is presented in Section El 



2 Post-Newtonian expansion 



We adopt the definition of a post-Newtonian approximation from , see also |22] for the 
Einstein case. Therefore matter and fields are described by functions (/^(c), i?(c), -B(c)) 
depending on a parameter c G [cq, oo) giving a one-parameter family of solutions of 
JRVMcD . This means that (/^(c),E(c),5(c)) describes a one-parameter family of solu- 
tions of physical systems which are represented in parameter-dependent units where the 
numerical value of the speed of light is given by c. A more conventional physical de- 
scription of the post-Newtonian expansion would say that in a fixed system of units the 
occurring velocities are small compared to the speed of light. To be more precise, note 
that (/^, E, B) is a solution of JRVMJ with c = e-^l^ if and only if 



B'{t,x) 



e^E{e^/\ex), 
e^Bie^^H^ex) 



(2.3) 



is a solution of dRVM, 
unchanged 



with c = 1. In this scaling the masses of the system remain 
J J f^'^dpdx = jj f^dpdx 



while the average momenta 



pf^'^ dpdx 



^11 



pf^ dp dx = y/ep 



are scaled by ^/e. By definition of the rescaled fields these fields are slowly varying in 
their space and time variables. Thus, the limit c ^ oo corresponds to an adiabatic limit. 
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In this work we treat the speed of hght c as a parameter and study the behavior of 
the system as c — * oo, but note that all Theorems can also be formulated in a parameter 
independent fashion. In that case the value of c is fixed, say c = 1, and the initial data 
has to be scaled according to ()2.3|) . see |5| for details. 

We start with a formal expansion of all coefficients occurring in ( RVMc ) in powers of 
c-\ 



E 
B 

P 
3 



fo + c-'ft + c-^f^ + c--'ft + ... 
Eq + c~^Ei + c~^E2 + c-^Es + ..., 
Bo + c^'Bi + c'^B2 + c~^Bs + . . . , 

PO + C^Vl + C~V2 + C"^P3 + • • • , 
Jo + C~^Jl + C~^J2 + C-^JS + ■■■■ 



(2.4) 



Moreover, also the initial denities are assumed to be expandable /°'^ = /g'^+c~^/{''^ + ... . 
While considering the Cauchy problem the same is understood for the electromagnetic 
initial fields. Finally p = p — (c~^/2)p^p + ... by where = These expansions 

can be substituted into ( RVMcD . Comparing coefficients at every order gives a hierarchy of 
equations for these coefficients. The equations in order k will be addressed as the fc/2-PN 



equations and the solutions (/±''=/2-p^ = Y/j=oC~^ ff ^ E^'"^ 



PN 



Qk/2 



PN 



.) as the 



fc/2-PN approximation contributing to the fact that in the context of general relativity 
post-Newtonian approximations are usually counted in orders of c~^. In order zero the 
well known Vlasov-Poisson system of plasma physics appears. 



5i/o^+P-V./o^±i?o-Vp/o^ = 



Eo(t,x) 



Po 



/o^(0,x,p) 



— j \z\ '^z po{t,x + z) dz, 

ifo - fo)dp, 
fo'^{x,p), 



(VP) 



where 



z\z\ 



Note that the degrees of freedom of the electromagnetic fields up to this order are lost, 
reflecting the fact that the limit c — > oo is singular and the hyperbolic fleld equations 
become elliptic. We recall that in it has been shown that as c — >■ cxd the solutions 
of JRVMeD approach a solution of the Vlasov-Poisson system with the rate 0(c~^); see 
[2 E] for similar results and JH] for the case of two spatial dimensions. The respective 
Newtonian limits of other related systems are derived in |28| IHj. 

Concerning a general k we assume that the lower order coefficients are already com- 
puted. Then the fields in order k have to solve 



VxEk 
Vx Bk 



—dtBk^i, 
dtEk-i + 4:7ijk-i, 



V-Ek 
V-Bk 



0. 



The corresponding Vlasov equation is given by 
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where Rk is depending on f^, '^xff, "^pf^ ^j) Bj, j = 0, • • • ,k — 1. If we assume that 
fl' — for all odd what we will do henceforward without mentioning, and employing 
the explicit form of Rk it is easy to show that we consistently can set 

/2±+, = 0, En+i^Q, B^i^O (2.5) 

for all / = 0, 1, 2, ... . Granted is known, we can easily solve for using characteristics. 
Note that for all orders k the characteristic flow is given by the vectorfield (p, ±£'o). Thus, 
for A; > 1 Ek only enters the Vlasov equation via the right hand side. 

Using the vector identity —V x V x +VV- = A we can rewrite the field equations. 

E2k = A7rA-\Vp2k + dtj2k-2) + /^~\d^E2k-2) 

k 

= 47r ^ [d? {Vp2k-2l + dtj2k-2l-2) 

1=0 

B2k+i = A-\d^B2k-2)-^7rA-\Vxj2k) 

k 

= -AnJ2^~'~'[df^xhk-2i], (2.6) 

1=0 

where quantities carrying a negative index are set to zero. For given f^, thus pk and 
jk are given, and if wc assume that all densities have compact support we can solve 
these equations using the convolution with the fundamental solution of A^^^'. Of course, 
without boundary conditions these solutions are not unique, and at least for higher orders 
these solutions will not vanish at infinity. 

Nevertheless, if we choose these fields the coupled equtions are easily solved by a 
fix-point iteration for E^. Thus, this (naive) PN approximation scheme is well defined. 

According to this scheme Bi is given by 

Bi{t,x)^ J \z\~'^z X jo{t,x + z)dz (2.7) 

where 

Jo = jp{U-fo)dp. (2.8) 
The couple (/2, £^2) is the solution of 

dtft +P- ^Ji -\P^P- ^Jo ± Eo ■ Vpft ±iE2+px B,) . Vp/o± = 0, 

E2{t, ~ 2 y ^^tP^^^^ X + z) dz — J \z\~^dtjo{t, x + z) dz — J \z\~'^zp2{t, x + z) dz 

P-^^ J if 2 - /2~) dp 

f2{0,X,p) = f2^{x,p). 

(LVP) 

In analogy to the particle model the 1-PN approximation 

f±,iPN^f±^^-2f±^ E^P^ = E, + c-'E2, B'^^ = c-'B, (2.9) 

is also called the Darwin approximation. It is Hamiltonian in the following sense. If the 
conserved energy 

S = Jj V'^ + c-V/2 (/+ + /-) dpdx+^j {E^ + S^) dx 



of JRVMel ) is expanded according to ()2.4|1 one obtains the Darwin energy defined by 



£d = £D,km + ^D,pot- The kinetic and potential energy are given by 

£D,kin = JJipV2-c-Ymf^ + fo)+c-y/2if+ + f^)dpdx and 

SD,pot = ^ I El + 2c-^Eo ■ E2 + c-^Bl dx 
ovr J 

respectively. One can check that £0 is conserved along solutions of the 1-PN approxima- 



tion. If in the Cauchy problem of ( RVMg ) we adapt the initial data ()1.2|) to suit the data 
of the IPN approximation the solutions are tracked down with an error of order c~^, see 
[5]. Hence, the naive post-Newtonian expansion is valid up to this order. 

2.1 Radiation damping in the 1.5 PN approximation 

Using the naive expansion, according to ()2.5p and (j2.6|) we would simply have to add 



53(t,x) = -47r A-i(Vxj2)+A-2(52Vxjo) H 



,X] 



2 

together with the factor c~'^ to the magnetic field, where 



z\'^V ^ ]2{t,x)dz+]- / \z\d'lV ^ ]Q{t,x)dz (2.10) 



32 = I pUt - f2) - jPif^ - fo) dp. 



Therefore the relevant energy So has not to be changed in comparison to the IPN order 
and we would remain with a Hamiltonian system. On the other hand it is known that 
in the full relativistic system energy is radiated to null infinity. In [HI Theorem 1.4] it is 
shown that in the limit c 00, corresponding to small velocities, the total amount of 
radiated energy is given by 

Ic-W, (2.11) 
where D is the dipole moment of the Newtonian limit of the matter defined by 

D{t) = jxpo{t,x)dx. (2.12) 

This theorem gives a mathematical formulation and a rigorous proof of the Larmor formula 
in case of Vlasov matter. Hence, we should introduce a radiation reaction force causing 
this loss of energy. As already suggested in ^3 HHI we modify the Vlasov equation of the 
Newtonian distribution by incorporating a small correction into the force term, 

dtf^ + P ■ ^Jo ± (^0 + ^D) ■ Vp/o^ = 0. (2.13) 

The additional term is the generalization of the radiation reaction force used in particle 
models, see jT^ formula (16.8)]. In passing we note that for this system the "energy" 

£s = l JJp\U + U)it,x,p)dpdx+^ J \Eo{t,x)\'dx-'^c-'D-D (2.14) 
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is decreasing, more precisely one obtains 



|£. = -|5|flWP, (2.15) 



the subscript 5* referring to the name "Schott" -energy under which this energy can be 
found in the hterature. This decreasing of energy can be attributed to the effect of 
radiation damping. If we remain with the positiv definite energy of the Vlasov-Poisson 
system 

£vp = \ Jjp\f^ + U)dpdx+^ j \E^{t,x)\^dx (2.16) 

evaluated along solutions of ()2.13|1 the "friction" has a definite sign only in the time 
average 

2.1.1 "Unphysical" solutions and the "Reduced Radiating Vlasov-Poisson 
system" 

Introducing system ()2.13|) one immediately runs into the problem that an initial datum 
has to be supplied for D (note that -D(O) and -D(O) are already determined by /q^) and 
there is no obvious way to extract this information from the approximation scheme. This 
phenomenon is also known in the theory of accelerated, and thus radiating, single charges 
and leads to the unphysical so-called run-away solutions. In JH] it has been observed that 
in the particle model this problem has the structure of a singular geometric perturbation 
problem, and the "physical" dynamics is obtained on a center-like manifold of the full 
dynamics. In order to adopt this language to the model under consideration here we 
assume that we are supplied with a (local in time) classical solution {f^,Eo) of ()2.13j) 
and assume that the support of f^(t,-,-) remains compact for all t in the interval of 
existence of the solution. We define the bare mass of the charges by 



M = JJ{f^ + fo){t,x,p)dpdx. 

Mass conservation and charge conservation for both species easily follow from ()2.13|1 and 
integration by parts, 

dtM = 0, dtp^ + \/.j^ = o^ 

where = J pf^ dp and Pg ~ I fo dp. We denote the additional degrees of freedom by 
y := D and compute with 

r] := (2/3)c"^M (2.17) 
and exploiting ()2.13j) in combination with integration by parts twice 

y = L)[21 + riD 

where D'^' is defined by 

D^^\t) = jj Eo{t,x){f^ + U){t,x,p)dpdx. (2.18) 
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Thus, we may rewrite ()2.13|) in a form clearly showing the structure of a singular pertur- 
bation problem. 

where and G are defined by 



(SGPP, 



Gif^,y)=y-D^'l 

In contrast to ^Hl we are dealing with a phase-space of infinite dimension. Thus, the 
proof of the existence of invariant manifolds is hard. We shall return to that question in 
a forthcoming paper. For the moment we shall take the existence of a smooth invariant 
manifold for granted and assume that it is given by means of a smooth function /i^ = 
h^ifS), acting on C^{R^ x R^) x C^(R^ x M^) and taking values in R^. In this subsection 
/q denotes the couple (/o'"*", /o'~)- The same convention will be used for /q = {/q, /g") and 
F = F-). The manifold Mrj = {(/q , /^^(/o )} is invariant under the flow of ( |SGPPJ if 



the solution of ( [SGPP,^ ) subject to the initial conditions (/o,i/(0)) = (/q, ^j?(/o)) satisfies 

y{t) = hr,ifo{t,;-)) 

for all times the solution exits. 

We want to establish a dynamics of Vlasov-Poisson type which is a good approximation 
of the dynamics on the manifold. For this reason we assume that we can expand /i^ in 
T] about 0, h^i = Hq + r]hi + 0{rf). Inserting 77 = we have G(/o, /io(/o! 0)) = which 
implies 

/io = /^['l. (2.19) 
Exploiting this information we find by a formal calculation 

^y = G{hMU)+rihi{h) + 0{r,^)) 
= vdyG{fo,ho{fo))h{fo) + 0{v') 

and on the other hand 

Vy = V<K{fo)Jo>+0{r]') 

which yields 

hifo) = [dyG{fo,hoifo))y' < K,F{fo,hoifo)) > . 
Here ' denotes the Frechet derivative and < ■, ■ > is the duality pairing, 

<Kifo),Fif,,hoifo)) >= Kifo)-jj {-p-^Mo~fo)-E,-VM + U)y;x,p)dpdx. 

In view of the definition of G and ()2.19|) we have 

y = Z)M+r/<DPl>(/o,Z^[^l)>+0(r^^). 
After some straightforward computations we find 
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where 



and 



Z}[=^](t) = 2 J H+{t,x)j^{t,x)-H-{t,x)j^{t,x)dx 
H^(t,x) := ^ \z\~^ { - 3z IS) z + id) p^{t, X + z) dz E W 



Note that 



(2.20) 
(2.21) 

(2.22) 
= 0. 

Ahernatively D'^' may be found by means of the following formal calculations, using 



Hiz) 



-3z ® z + id 



is bounded on \ {0} homogeneous of degree zero and satisfies J\^\^i H{z) ds{ 



j Ut,x)dx + 0{c-^) 
= Z^'Pl + C(c-3) 

Whereas this second derivation is more simple the first derivation reveals the import 
connection to the dynamics on the manifold. 

We introduce the "reduced radiating Vlasov-Poisson system" 

dtf^ + P ■ VJ^ ± {Eo + 2/3c-3^[3l) ■ Vp/o^ = 0, 
Eo{t,x) = — J \z\~'^z pQ{t, X + z) dz, 



Po 

/o^(0,a;,p) = /o'^(x,p), 



(rrVPe) 



where D'^' is defined according to ()2.2()j) and 1)2.211) . 

Our first proposition addresses the existence and uniqueness of local classical solutions 
of ( rrVPc ). Furthermore, it provides us with some estimates useful in the sequel of this 
paper. In every (even) order k of we shall only consider smooth compactly supported 
initial data, 



(2.23) 



fl'^{x,p) = for |a;| > tq or \p\ > tq \\f^^\\w^,^ < So 



with some tq > and 5*0 fixed. We call the constants tq and 5*0 'basic' as all bounds 
related to {f^,Eo) occurring in the sequel will only depend on and Sq. The following 
proposition is proved in 3 . 

Proposition 2.1 // /q'^ satisfies \2.2'J^) then there exists a constant < T < oo such 
that the following holds for c > 1. 

(a) There is a unique classical, i.e. solution (f^jE^) of HrrVPc ) existing on a time 
interval [0,Tc) with T <Tc < oo. 

(b ) For every T < T there is a constant Mi (T) > such that for all < t < T 

f^{t,x,p) = if \x\ > Mi{T) or \p\ > Mi{T). 
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(c) In fact, is C°° and for every T < T there is constant M2{T), such that for all 
0<t<T 

\d'^f^{t,x,p)\ + \d^Eo{t,x)\ + \d^D^'\t)\ < M^iT) 
for all X eM?, p e M^, |a| < 4 and (3 <2 and 7 < 1 

(d) For this solution {fQ,Eo) we have 

where -D^^l is defined according to \2.1^) . 

Note that the constants T , Mi{T) , M2(T) appearing in Proposition 12.11 do only depend 
on ro and So, in particular they are independent of c. 

As the second moment JJ p'^ifo + fQ)dpdx cannot be bounded a priori by using 
energy conservation it seems difficult to prove global existence of classical solutions of 
( rrVPcD . Note that both methods yielding global existence of Vlasov-Poisson type systems 



essentially rely on such an a priori bound, see [21] or [23] and |20j. 

We shall use solutions of (rrVPc) instead of solutions of (jVPj) in order zero of our 



post-Newtonian approximation and define Bi, (/2 , E2) and B3 according to ()2.7|1 . (jLVPj) 



and ()2.10|) respectively. It is important to note that (/o, Eq) is the solution of ( rr VP^ ) . In 
particular this solution is depending on c and thus all other quantities defined here are 
depending on c. Concerning the solvability of ()LVP|) we have the following lemma whose 
proof is sketched in [Hj. 

Lemma 2.2 Let fl'^, k = 0,2 satisfy i2.2!^} . Suppose that {fQ,EQ) is the solution of 
j^rrVPc]) . Compute Bi according to \2. Then \LVF]) has a unique classical solution 
(f^,E2) existing on [0,Tc) and enjoying the following properties. 

(a) For every T < T there is a constant M^{T) > 0, such that for all < t < T 
f^{t,x,p) = if \x\>M^{T) or |p|>M3(T). 

(h) In fact, f2 is C°° and there is a constant M^{T), such that for all < t < T 
\d''f^{t,x,p)\ < Mi{T) for all xeM^pGM^ and \a\ < 2 

In the following section we want to show that 

E^ = Eo + c-^E2 + 2/3c-'^D^^^ (2.24) 
= c-^Bi + c-^B3 



yields a higher order pointwise approximation of (j^VMJ) than the Vlasov-Poisson or 
the Darwin system defined in [^j. We call ()2.24j) the radiation approximation. In the 
terminology of post-Newtonian approximations it is the 1.5PN approximation. Employing 
the Vlasov equation and integration by parts it is not difficult to prove the following 
formulas. 
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Corollary 2.3 The fields and can he written as 
E^{t,x) = - J \z\^'^z{po + c~'^P2)it,x + z) dz 

+c"^^ jj \z\'^!^3{z ■ pfz - p^z'^iftl' - fQ)(J:,x + z,p) dpdz 

jj \z\~'^ ® z + 1 

[Eo{t,x + z) + c-32/3DP](t))(/+ + /o-)(t,a; + z,p) dpdz 

+c-=^-D[=^](t) (2.25a) 
3 

and 

B^{t,x) = c-^ jJ \z\-h Ap{f^ - f^){t,x + z,p)dpdz (2.25b) 
-c"^^ jj \z\-^iz-pYzApif+-fo)it,x + z,p)dpdz 



+ ^ JJ \z\''^^{z Ap)(^z + {z-p)z A{---) 

(Eo(t, x + z)+ c~'2/3D^'\t)') (/o+ + fo){t, x + z,p) dpdz 

-c'^ j zA {H^{t,x + z)jQ{t,x + z) - H^{t,x + z)j^{t,x + z)) dz. 

Closing this section we want to mention that there is another variant of a damped Vlasov 
Poisson type system investigated in ^H] and ^E] ■ While for that system a global solution 
theory is at hand, the author did not succeed in comparing approximations based on 
solutions of that system on the one hand and solutions of the full system on the other 
hand. 



3 Comparison of the 1.5 PN dynamics with the Vlasov- 
Maxwell dynamics 

3.1 The Vlasov-Maxwell dynamics with adapted initial values 



For achieving the improved approximation property we match the initial data of dRVMc ) 
by the data for radiation system. For prescribed initial densities /^'^, = 0,2 we are 
able to calculate {f^,Eo), Bi, (/^,i?2) and B^ according to what has been outlined in 
Section |21 We then consider the Cauchy problem of ( RVMe ) where the initial values are 
given by 

fHO,X,p) = f°'Hx,p) = /o'^(x,p) + C-'f°'^{x,p) + C-V:jtee, ] 
E{0,X) = E°{X) ■= Eo{0,x) + C-^E2{0,X) + |c-3D[31(0) + C-^Jree ( (1^) 

B{0,x) = B°{x) := c-'Bi{0,x) + 0-3^3(0, x) + c'^S"^,,,. J 

In contrast to the contributions in the orders to 3, which are fixed by the values of 
the approximations, (/°'^ee' -^cjree' ^cjree) be chosen freely only subject to the 
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constraints V ■ E°j^^^ = An - f°j^J dp and V ■ B°f^^^ = 0. Note that the 

constraint equations in the lower orders are satisfied by fiat. Furthermore, we shall assume 
that the following bounds hold uniformly in c. 

f°cfree = ^ if 1^1 >^0 Or |p| > Tq (3.26) 



' c,free 

' c,free I 

\EcJree\\w^'°° + \\Bc,free\\w^-°° ^ 'S'o 



\\fc'tree\\L°° < 'S'o 



Before we formulate our theorem, let us recall that solutions of (|RVMe| ) with initial data 
flK^jl exist at least on some time interval [0,T), which is independent of c > 1; see fl^ 
Thm. 1]. 

Proposition 3.1 Assume that fl'^.k = 0,2 satisfies KTM) . If f°^^ , E° and B° are 

defined according to 1/6]) . then there exits < T < oo ( independent of c) such that for all 
c > 1 the system \RVMc ) with initial data 1/6]) has a unique -solution (/, E, B) on the 



time interval [0,T). In addition, for every < T < T there are constants M^{T)^ Mq{T) 
(independent of c) such that 

f^{t,x,p) = if \x\>M,{T) or \p\>M,iT), (3.27a) 
\f^{t,x,p)\ + \Eit,x)\ + \Bit,x)\ < M6(T), (3.27b) 

for allx, p eR^, te [0,T], and c > 1. 

Actually, in |24j Thm. 1] E° and B° do not depend on c but an inspection of the 
proof shows that the assertions remain valid for initial fields defined by After these 
preparations we can state the first of our main results. 

Theorem 3.2 Assume that f^'^,k = 0,2 satisfies \2.2'J^) . From fl'^ calculate {fQ.Eo), 
Bi, (/^,£^2) CLnd B^, choose initial data {f°'^, E°, B°) for HRVMcD according to 1/6]) 
and llS.2f^} . Let (/, E, B) denote the solution of \RVMc[ ) with initial data 1/6]) and let 
{f^,E^,B^) be defined as in \2.2J^ . Then for every T < min{T, T} and R > there are 
constants M{T) > and M{T, K) > 0, such that 

\f^it,x,p)-fi{t,x,p)\ < M{T)c-' (xGM^), 

\E{t,x) ~ E^{t,x)\ < M{T,R)c-^ {\x\<R), (3.28) 
\B{t,x) - B^{t,x)\ < M{T,R)c-^ {\x\<R), 

for all p G M^ t G [0, T], andc>l. 

The constants M(T) and M(T, R) are independent of c > 1, but do depend on the basic 
constants rg, S'o. Note that if dRVMeD is compared to the Vlasov-Poisson system ()VP|) 
only, one obtains the estimate 

\fit,x,p) - fo{t,x,p)\ + \E{t,x) - Eo{t,x)\ + \B{t,x)\ < M{T)c-\ 
see PH, Thm. 2B], and if compared to the Darwin system the estimates 

\f{t,x,p) - r(t,x,p)| + \B{t,x) - i?^(t,x)| < M(r)c-3 



13 



\E{t,x) - E^{t,x)\ < M{T,R)c-\ 

see 13 Thm. 1.1]. On first glance it could seem that the time interval [0, T] C [0, min{T, T}), 
which might be very small, is a strong limitation of the theorem. But if the theorem is 
formulated in an ^-depending fashion using fixed units as indicated in Section |21 and elab- 
orated in 1^1 the approximation is valid on the time interval [0,e~^/^T] and thus for long 
times on that time scale. 



3.2 The retardet Vlasov-Maxwell dynamics 

Following P we introduce the retardet relativistic Vlasov-Maxwell system. 

dtf^ + P ■ V./± ±{E + c-'p X B) ■ Vp/± = 0, 



E{t,x) = 
B{t,x) 
P = 



dy 



c 



x-y 

-1 



(Vp + c '^dtj){t-c ^\x-y\,y) 
dy 



\x - y\ 



-Wxj{t-c \x-y\,y) 



if^-Ddp, 



P{f^ - f )dp 



RVMJ 



If we assume that is a global solution of drctRVMj and that also E and B are C^, 
then, by means of the Vlasov equation, p and j satisfy the continuity equation 

5iP + j = 0, (3.29) 

and therefore the retardet fields are a solution of Maxwell's equations. Thus, (/, E, B) 
also solves ( RVMe ) . Note that it is necessary to know the densities for all time (— cxo, t] in 
order to compute the fields at time t. Hence, there is no sense in the notation of a local 
solution of this system. As in the case of the Cauchy problem every solution of jretRVMc ) 
satisfies the identity 



= /±(0,X±(0;t,x,p),P±(0;t,x,p)), 
where s i— > (X=''(s; t, P^(s; t, solves the characteristic system 

X = p, p = ±{E + c^^p X B), 
with data t, x,p) = x and t, x,p) = p. For this reason 

o</^(t,x,p)<||/±(o,-,-)IL. 



(3.30) 



(3.31) 



(3.32) 



In order to derive our results on PN-approximations let us assume that we are furnished 
with a one-parameter familiy of solutions (retRVMc) satisfying certain plausible a priori 
bounds and smoothness conditions. We consider solutions of jretRVMj passing through 
a certain density configuration at time t = 



c,free 



(3.33) 



where /^'^, k = 0,2 and /°'^ee satisfy (j2.23|l and (j3.2(i|l respectively with some constants 
To and Sq. 
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Assumption 3.3 (a) For each c > 1 there is a solution e C"^(M x x M^) of 
HretRVMc ), passing through /°'^ at time t = 0, i.e. 



(h) There exists a constant Pi > such that f'^{t,x,p) = for \p\ > Pi and c > 1. In 
particular, f'^(t,x,p) = for \x\ > Tq + by \3. 31]) . 

(c) For every T > 0, R > 0, and P > there is a constant M-j{T, R,P)>0 such that 

\d:+'f^{t,x,p)\ + \d'^VJ^{t,x,p)\ < Mr{T,R,P) 

for \t\ < T , \x\ < R, \p\ < P, and a = 0, 3, uniformly in c> 1. 

After these preparations we can state our second main result. 

Theorem 3.4 (Approximation of the retardet Vlasov-Maxwell system) 

Assume that f°'^, k = 0,2 satisfies WJ^) . From f°'^ calculate {f^, Eq), Bi, (f^, E2), 
B3 andD^^l D efine (/^, E^, B^) according to A ssume that (/^, i?, B) is a family 

of solutions of ^^^^M^ satisfying Assumption \3. 'A with constants Pi and Mi{T, R, P) . 
Take T from Proposition \2.1\ Then for every T < T and R> there are constants M{T) 
and M{T, R) such that 

\f^it,x,p)-fi{t,x,p)\ < M(T)c-^ (xGM^), 

\E{t,x) - E^{t,x)\ < M{T,R)c-^ {\x\<R), (3.34) 
\B{t,x) - B^{t,x)\ < M{T,R)c-* (k|<P), 

for allpe M^ t G [0,T], and c>2Pi. 

The constants M{T) and M(T,R) do only depend on rQ,So,Pi and Mi{-,-,-). In 
particular they are independent of c> Pi. 

3.3 Comparison with the particle model 

We shall compare our results for the Vlasov model with the results for the particle model 
governed by the Abraham system obtained in ^\ and ^Hj- Both systems share the 
features of Hamiltonian approximations up to 1-PN order and dissipative corrections in 
the 1.5-PN approximation leading to an increase of the phase-space. The right comparison 
dynamics is given on a center-like manifold, where the dynamics on this manifold can 
be approximated by a modified Vlasov-Poisson equation and a second order equation, 
respectively. In 16, section 3] it is shown that the force terms of the 1.5PN approximation 
used here do agree with the infinte particle limit of the comparison dynamics used in jTH] • 
In contrast to the PN approximation in this paper in ^\ and [THj only the forces 
are expanded but the main difference is in the treatment of the initial data. For the full 
particle model the initial data for the fields are supposed to be of "charged soliton" type. 
One can think of these fields as generated by charges forced to move freely for — 00 < t < 
with their initial velocity. For the approximation this leads to an initial time slip to which 
the charges need to "forget" their initial data. The initial data of the approximation is 
fixed by matching the data of the full system at time to- Therefore the initial data for 
the approximation are given only implicitly, first one has to compute a solution of the full 
system over a time span tg- Regarding the Cauchy problem of Vlasov-Maxwell system 
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we do the matching the other way round. For a given initial density one computes the 
fields of the approximations and imposes their values at t = as initial data on the fields 
of the full system. Therefore these initial data are given more explictly. Even more it is 
possible to calculate them by the values of /q'^ and /a'^ alone, see ()2.25a|) and ()2.25b|) . 

Moreover both results in this paper (Theorem 13.21 and Theorem 13. 4|) seem to be 
stronger than the results obtained for the particle model. In ^Hj the error bounds of 
the 1.5-PN approximation in contrast to the error bounds of the 1-PN-approximation are 
only improved in a certain direction, see [HI formulas (3.21) and (3.32)]. It seems rea- 
sonable that a matching of the initial conditions at time t = according to the treatment 
of the initial conditions used here might improve those bounds. 



3.4 Notation 

In the remaining more technical sections of this paper we shall make use of the following 
notations. -8(0, R) denotes the closed ball in with center at x = or p = and radius 
R> 0. We write 

g{x,v,t,c) = Ocptic-""), 

if for all i? > and T > there is a constant M/j > only depending on the basic 
constants ro,So and, while dealing with solutions of jrctRVMcT ), Pi and Mt{-,-,-) such 
that 

\gix,v,t,c)\<Mnc-"' (3.35) 
for |x| < R, p eR^, t e [0,T], and c > 1. Similarly, we write 

g{x,v,t,c) = 0{c-n, 

if for all T > there is a constant M > only depending on the basic constants such 
that ()3.35p holds for all x,p E M^, t G [0, T], and c > 1. In short, O is global in x,p and 
c while Ocpt is local in x global in P and c, and both symbols are local in t. In general, 
generic constants only depending on the basic constants are denoted by M. Furthermore, 
in the following sections we shall use the notation 

f = r-f-, fi = r,^-fi''. = /rand fi = f^-fr 

for k = 0,2 and / = 0,2, R. (Recall that /g and /o were used in a different way in 
Subsection 12.1.11 ) 



4 Proof of Theorem 13.2 
4.1 Estimating E — 

In Section 16.11 below we will show that the approximate electric field from ()2.24j) 
admits the following representation. 

= El + + + Coptic-') (4.36a) 

with 

E^^,{t,x) = [ { -\z\-^z{po + c-^P2)-c-^\z\-'dtjo (4.36b) 

J\z\>ct L 
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+ l/2c-'z9iVo + [^o(Po + Po )] a; + z) dz 

El^{t,x) = I [ \z\-^K,iz,c~'p)fRiiiz),x + z,p)dpdz (4.36c) 

J\z\<ct J 

/ / \z\-^K2iz,c-^p)Eo{i{z),x + z)Ur + fR){i{z),x + z,p) dpdz 

J\z\<ct J 

ELndit,x) = (ct)-' [ [ K,{z,c-'p){f° + c-'f°)ix + z,p)dpds{z) (4.36d) 

J\z\=ct J 



— l/3tc / zz ■ d^jo{0,x + z) ds{z) 

J\z\=ct 

+2/3c~^ / {z ■ p)dtjo{0,x + z,p) dpds{i 

J \z\=ct J 



where the subscripts 'ext', 'int' and 'bound' refer to the exterior, interior and boundary 
integration in z. The kernels are given by 

Ki{z,p) = —z + 2z{z-p)—p — ?)z{z-pY + 2{z-p)p + zp^ + Az{z-pf 

-3p{z -pY - 3z{z ■ p)p^ + 3/2pp2 (4.37a) 
K2{z,p) = z ® z — 1 — 2z ® z{z ■ p) + z ■ p + p ® z (4.37b) 
K^{z^p) = z{z ■ p) + p{z ■ p) — z{z ■ pY + z{z ■ pY — p{z ■ pY — l/2z{z ■ p)p^ . 

We also recall that z = z\z\~^ and i(z) = t — c~^\z\. On the other hand, according to 
Section Ell below, we have 



E = E,,t + Eint + E,ound + 0(c-') (4.38a) 

with 



Ee.t{t,x) = / I -\z\-'z {l + tdt){po + c-'p2) + (l/2t'dt + l/Qrdt)p^ {A.?>^h) 



\z\>ct 

+ll2c-^z{dl + t9f)po - c~M'\dt + tdt)jo 
+2/3c-39t [Eo(p+ + Po )] } (0, x + z)dz 
Eint{t,x) = / \z\'''^Ki{z,c~'^p)f{i{z),x + z,p)dpdz (4.38c) 

J\z\<ct J 

\z\~^K2{z, c-^p)E{i{z),x + z){f+ + f-){t{z),x + z,p) dpdz 



z\<ct 



Ebound{t,x) = E^^Ut.x)- (4.38d) 

In order to verify ()3.28p . we fix constants i? > and < T < min{T, T}. For 
X G -8(0, R) and t < T we start by comparing the exterior fields. We obtain from ()4.38b|) 
and ()4.36b|) . due to \z\ = 1, and taking into account Theorem I2.1f b)-fc) and Lemma 
I22i;a)-(b) 

\E^,,{t,x)-E,,tit,x)\ 

< [ \\z\-^\poit,x + z)-{l+tdt + l/2t^d^ + l/Qt'^df)poiO,x + z)\ 

J\z\>ct I 
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+c~'^\z\-^\p2{t,x + z) - {l + tdt)p2{0,x + z)\ 
+c-^\z\-' j \p\ \dth{t,x + z,p) - (1 + tdt)dtfo{0, x + z,p)\dp 
+1/2 c"ViVo(t, x + z)-{l+ tdt)d^po{0, x + z)\ 
+2/3c-3|9i[Eo(p+ + Po )](t,x + z) - dt[Eoip^ + Po )](0,x + z) \ dz 



< MM2{T)Ml{T) [ \z\-HHBio,M.{T))ix + z)dz 

J\z\>ct 

+c-^MM^{T)MliT) [ \z\-hHBio,AMT)){x + z) dz 

J\z\>ct 

+c-^MM2iT)Mf{T) [ \z\'HHB(^o,AMT))ix + z)dz 

J\z\>ct 

+c-^MM2{T)Ml{T) [ t^lBio,AMT)){x + z) dz 

J\z\>ct 

+c-^MM^{T)M^iT) [ tlBio,MiiT)){x + z) dz 

J\z\>ct 

< m[ {t^\z\-^ +t^C-Wz\~^ + \z\~' + l)+tc'^)lBiO,R+Mo){z)dz 

J\z\>ct 

< Mrc-^- (4.39) 
where Mq = m.a2c{Mi{T) , M^{T) , M^{T)} and we used that for instance 

[ \z\'HBio,B+Mo)ix + Z)dz< {ctyH" [ \z\^dz < MrC-\ 

J\z\>ct J\z\<R+Mo 

To bound \E\nt{t^x) — E^^^{t^x)\^ we first claim that 

\E{t,x)-E^{t,x)\ = 0{c-^). (4.40) 

which can be proved analogously to Theorem 1.1]. (This estimate holds uniformly in 
X G M^.) Next we define 

g±(t) = sup{|/±(s,x,p) - /|(s,x,p)| : X G M^p G M^ s G [0,t]}. 

and Q{t) = Q^{t) + Q^{t). Recall that f^{s^x,p) = f^{s,x,p) = /^(s,x,p) = if 
> Mq or \p\ > Mq. Therefore the kernels Ki, i + 1, 2, 3 are bounded on the domain of 
the densities by M(l + M§) uniformly in c > 1. From ()4.38d;i . (j4.3(icjl . (jOOl), (|E3), and 
< t{z) < t for \z\ < ct we obtain 

\E,^t{t,x)-El^,{t,x)\ 

< M(i + Mo)' I j [{W + c'M-'\E,\) {\f^ - ft\ + ir -/«!)) 

J\z\<ct J 

+c-M~^\E - Eo|(|/+| + ir I)] {i{z),x + z,p) dpdz 
r3\/\/r3/i , T,jtrr\\rM+\ I i Ul-l 



< M(l + Mo^)Mo^(l + M2(T))Q(t) / {\z\-' ^\z\-^)lBio.M,){x ^ z)dz 

J\z\<ct 

+c-^M(l + Ml)MlM^{T) [ \z\-' lBio,Mo)ix + z) dz 

J\z\<ct 
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< MR{c-^ + Q{t)), (4.41) 
since for instance 

/ \z\~^1b{o,Mo){x + z)dz < / Izl'"^ dz < Mr. 

J\z\<ct J\z\<R+Mo 

Reca lling that the Ebound{t,x) = E^^^^^{t,x), we can summarize ()4.H8aj). ()4.aeiaj]. (jOHI), 
and (|OT|) as 

\E{t,x)-E'^{t,x)\ <Mn{c-^ + Q{t)), (4.42a) 

for \x\ < R and t G [0, T]. Using (jfi.fill) and (jdC^I) below and similar arguments as for the 
electric fields yield 

|5(t,x) -5^(t,x)| <MRic-^ + Qit)), (4.42b) 
for |x| < R and t G [0,T]. 

4.2 Estimating /± - /| 

It remains to estimate = - f^. Using JRVMcD , (jT^ . JrrVPJ , and (fnVP| . it is 
found that 

+ p ■ V. ± (E + c-^p X 5) • Vp 
= -dJi-p-VJ^T{E + c~'pxB)-Vpf^ 
= (p- 1/2 c- V p - p) ■ V./o^ + c-2(p - p) . VJt 
±{{E^ -E) + c-'p X (5^ - B)) ■ Vp/I 
±c-i(p - p) X 5^ ■ Vp/o^ ± c-V X i?3 ■ Vp/o^ 
Tc-2(E« - Eo) • Vjt T c-'p X 5^ ■ Vp/2±. 

If IpI < Mo, then also |p| = (1 + c"'^p'^)~^^'^\p\ < \p\ < Mq uniformly in c, and hence 

\p- {l- l/2c~V)p| < Mc'^. 

Next we note the straightforward estimates \B^{t,x)\ < Mc^^, \E^{t,x) — EQ{t,x)\ < 
Mc~'^ and |i?3(t, a;)| < M for t G [0,T]. In view of the bounds in Theorem I2.ir c) and 
Lemma |22tc), thus by (|4.42a,jl and ()4.42bjl . 

\dtq^{t,x,p) + p -V ^q^{t,x,p) + {E{t,x) + c-^p x 5(t, x)) ■ Vp g^(t, x, p) | 

< Mmo(c-^ + Q{t)) = M{c~^ + Q{t)) (4.43) 

for |x| < Mo, \p\ < Mo, and t G [0,T]. But in {{t,x,p) : |x| > Mo} U{(t,x,p) : \p\ > Mq} 
we have = — = hj the above definition of Mo > 0. Accordingly, ()4.43p is 
satisfied for all x G M^ p G ]R^ and t G [0, T]. Now, for any x G R^ p G and t G [0, T] 
we compute using ()3.3()j) and ()3.31|1 as well as ()4.43j) 

^g±(s,X±,P±)| = |9ig±+p-V,g±±(E + c-ipx5)-Vpg±||(,,x±,P±) <Mc-^ + MQ(s), 
for s G [0,r]. Here the characteristics are evaluated at {s;t,x,p). Note that 

g±(5,X±,P±)|.=o = (/^ - /^)(0,X±,P±)|.=o = C-^(/e:>:ee " /cT.ee) (^^ l.=0- 
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Thus, 

\q^{t,x,p)\ 



g±(s,(X±,P±)(s;t,x,p))||,=i 



ds 



< c-V^freeiO, (^^ ^^)(0; t, x,p))\ + | (X±, P±)(.; t, 

< Mc-^ + M [ Q{s)ds. 

Jo 

But by the definition of Q{t) it follows that 

Q{t) < Mc-^ + M [ Q{s)ds 
Jo 

for t G [0,T]. Then Gronwall's inequality implies Q{t) < Mc~^ for t G [0,T]. Inserting 
this into ()4.42a|l and (|4.42b|l yields the assertion of Theorem . 21 and completes the proof. 
□ 



5 Proof of Theorem 13.41 

We recall the following representation of the electric field E from jH]. 

E{t,x) = jj \z\-'^Kt{z,p) f{*,p)dpdz (5.44) 
jj \zr'Ks{z,p){E + c-'p A + dpdz 

where the integrals are to be extended over M.^, Kt and Ks are given by ()6.63a|) and 
fl6.63b|) . respectively and 

(*) = {t — c~"^|z|, X + z). 

Now we fix some constants T < T and R > 0. Furthermore, we assume that c > 2Pi as 
well cLS I X I < R and \t\ < T. Define 

p* = p*{T) = max{Pi, Mi(T), M^iT)} (5.45a) 

and 

r* = r*{T, R) = max{2(ro + PiT) + P, Mi(T), MsiT)} (5.45b) 
Then |x + 2;| > r* implies 

1 1 

To + Pi\t — c^"'"|-2| I = ro + Pi|t — c~^\x + z — x\ \ < ro + PiT + -R + - \x + z\ < \x + z\. 

Therefore, \p\ > p* or |x + 2;| > r* yields f^{*,p) = by Assumption Kl^f b). This 
argument shows, that in ()5.44|1 we can replace // dpdz by J^^^^^^j.* i|p|<p. dpdz. In other 
words, in these integrals we may always assume that both \z\ and \p\ are bounded by a 
bound depending on tq. Pi, R and T. Since the p-domain is bounded we have 

Kt{z,p) = Ki{z,c-^p) +0{c-^) (5.46a) 
Ks{z,p) = K2{z,c-^p)+0{c-^) (5.46b) 
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where Ki and K2 are defined in ()4.37ajl and (j4.37bjl . respectively. Furthermore, 

2 

c-^p = c-^p- c~^^p + 0{c-^). (5.46c) 
Firstly, we have using Assumption I3.3f b).fc) 



\B{t,x)\ = c"W |^|-^|V X (5.46d) 

J \x+z\<r* 

< c-'M7iT + R + r*,r*,p*) [ \z\~' dz [ dp = 0,pt{C^). 

^|a;+2j<r* \p\<P* 

In the same waywe are able to bound the electric field 



\d\E{t,x)\ = j \z\-^ d\S/p + c-^dl-^^] {*)dz (5.46e) 

J \x+z\<.r* 

< M-r{T + R + r*,r*,p*) / \z\-Uz / dp = Ocpt{l) 

J\x+z\<r* J\n\<n* 



'\x+z\<r* "J \p\^P* 

for / = 0, 3. Hence, using fl5.46d|) 

j j \z\-'Ksiz,p)c-'pABi*)ir + f-)i*,p)dpdz = 0,ptic~^). 
For this reason and using ()5.46e|) we conclude for the second term in ()5.44p 

JJ \z\'^Ks{z,p){E + c-'p A + n{*,p) dpdz 

= c-^ JJ \z\-'Ks{z,p)E{*){f+ + f-){*,p)dpdz + 0,pt{c-^) = 0,j,t{ 



Secondly we shall expand the retardet time. For every smooth function ip we have 



|,-y|2 |»/|3 |^|4 



with some t — c~^|2;| < C, < t. Using Assumption I3.3f b).(c) and ()5.4fie|l we have 

/ kr'+' / KT{z,p)dlf{^{t,x,z,p),x + z,p)dpdz = 0,pt{l) (5.46f) 
for / = 0, . . . , 4 as well as 

/ \z\-'^' I KL{z,p)dl[E{f^ + n]{at,x,z,p),x + z,p)dpdz = 0,pt{l) 

J\x+z\<r* J \p\^P* 

(5.4 6g) 

for / = 0, 1,2 and any choice of t — c~^\z\ < ^(t,x,z,p) < t Now from ()5.46a|) -( '5.46g| ) we 
can expand the electric field in powers of c~^. To order zero we have 



E{t,x) = - \z\~'^zp{t,x + z)dz + Ocptic~^) 

J \x+z\<r* 

= - J \z\''^zp{t,x + z) dz + Ocpt{c~^)- (5.47) 
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Here we employed that f^{t, x + z, p) = if \x + z\ > with 

rt = ro + PiT, (5.48) 

see Assumption 13. 3r b) and use r* > r\ see (j5.45b|l . In the first order we obtain two terms, 
the first coming from the expansion of the kernel Kt and the second from the expansion 
of the retardet time, 



and 



\z\ '^{2z{z ■ p) — p)f(t,x + z,p) dpdz 



Here and in the sequel we shall treat terms containig time derivatives in the following 
way. Firstly, we replace dtf"^ by —p ■ Vxf"^ T {E + c'^p A B) ■ Vp/^ using the Vlasov 
equation. Secondly we do an integration by parts utilizing 

{E + c^'p A B) ■ VJ^ = Vp ■ [{E + c-'p A B)f^] . 

For the second term we therefore obtain 

—c^^ // \z\^^{2z{z ■ p) — p)f{t,x + z,p) dpdz. 



Hence, using ()5.4(jc|l . the contribution in first order vanishes and in ()5.47|1 we can replace 
Ocpt{c~^) by (9cpt(c~^), but we have an additional term occuring in the third order, namely 

JJ \^\~'^p'^i^i^ ' P) ^ fi'^^^ ^^P) dpdz. (5.49) 

In the next orders we will simply replace p by p according to ()5.4(jc|l without comment. 
For the second order we start with the terms with two time derivatives; 

-c"^iy" zd^p{t,x + z)dz (5.50a) 

= -c-^^ jj \z\-^{-z{z-p)+p)dtf{t,x + z,p)dpdz + 0{c-^). 

Now we collect and rewrite the terms with one time derivative, including the term coming 
from ()5.50a|) . 



\z\ ^ (^z{z ■ p) — p + -{—z{z ■ p) + p)^ dtf(t, X + z,p) dpdz 
/ 9 3 \ 

i — 2^(^ ■ P)"^ + ■ p)p+ 2^P^ j fi^i X + z,p) dpdz (5.50b) 

- jj \z\-^ (^^ ® ^ - ^(^' ^ + ^)(/^ + /")('^' ^ + z.p) dpdz + C(c-^). 
Together with the remaining terms coming from the expansion of we therefore obtain 
E{t,x) = jj \z\~'^i^- z + c~'^(^^z{z ■ pf -^zp^^^f{t,x + z,p)dpdz (5.50c) 



22 



-c-'i jj \z\-^ z + l^E{t,x + z){f+ + f-){t,x + z,p) dpdz + Ocptic^'^) 

At last we turn to the third order and, following the usual route, first treat the term with 
three time derivatives, 



c~^^ I zd^p{t,x + z)dz = c-^^ Jj pd'ff{t,x + z,p)dpdz + 0{c-^). (5.51a) 



Regarding terms containing second time derivatives including the term coming from 
fl5.51aj) we have 

' ' ^ z{z ■ p) p -\ — p\d1f{t, X + z^p) dpdz 

2 6 



= c ^ jj \z\ ^ <y— 2z{z ■ pY + [z ■ p)p + zp^ jdtfit^x + z,p) dpdz 

Now we collect those term which contain exactly one time derivative including those 
coming from (j5.51bj) . 



c jj \z\ ^ — {—3z{z ■ p)'^ + 2{z ■ p) p + z p"^) — 2z{z ■ p)'^ + {z ■ p) p + z p 
dtf(t, X + z,p) dpdz 
+^"^ jj l^-iz(^z-l) + z®z- ^^dt[E{f+ + f-)]{t,x + z,p)dpdz 

= c^^ // |-2|~^< — Az{z ■ pY + 'd>{z ■ pY p + 2{z ■ p) p^z — p^ p> f{t,x + z^p) dpdz 



+c ^ <2{z ■ p)z z — p z — z ■ p>E{t,x + z){f^ + f ){t,x + z,p) dpdz 



+^~^1^t jj E{t,x + z){f^ + f-){t,x + z,p)dpdz + 0{c-^). (5.51c) 

Note that we do not touch upon the time derivative in the term of the last line. This term 
is responsible for radiation effects. Collecting all terms without time derivative coming 
from the third order in the expansion of Kt, the first order of the expansion of Ks, 
()5.51c|) and ()5.49|) we note that all these terms cancel exactly. Thus, the only remaining 
contribution in third order is the radiation term 

dt jj Eit,x + z)if+ + f-)it,x + z,p)dpdz = dt jj E{t, z){f^ + f-)it, z,p) dpdz =: Bit) 

Next we try to recast this term making it comparable to D^^\ Note that using the 
bounds on see Assumption Kinf c). and the bounds on the support of with some 

t - c~'^\z\ < ^{t,x,z) < t 

r d7 

dtE{t,x) = -dt {Vp + c-^dtj){*) — 



\z\ 
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(a,vp + c-^95-)(*)n 
5*vp(*)^ + apt(c-^) 

|Z + 2|<T-* 1^1 

dz 



j dtVp{t,x + z)- 



\x+z 

-1 



+c-^ I dfVp{^{t,x,z),x + z)dz + Ocpt{c-^) 

\x+z\<r* 

= - J \z\''^zdtp{t, x + z)dz + Ocpt{c~^) 

Furthermore, we recall ()5.47p and emphasize that the continuity equation ()3.29|) holds for 
both species seperately. Using these ingredients we compute 

D{t) = dt [ [ E{t,z){f^ + f-){t,z,p)dpdz 



-ai(p+-p-)(t,y)(p+ + p-)(t, z) 



ip+-p-){t,y)dt{p+ + p-){t,z) 



dydz + 0{c-^) 



-2 // \y-z\-\y-z) 



dtp'^it,y)p {t,z)-dtp {t,y)p+{t,z) 



2 j {H+{t,z)r{t,z)-H-{t,z)j+{t,z))dz + 0{c-'), 



dydz + 0{c"^] 



where 



\z\~^{-3z 0Z + l)p^{t, x + z)dz 
and p^ and are defined in the obvious way. Summarizing, we have 



x) = 



— z + c '^(^-z{z ■ — -zp'^^jf{t,x + z,p) dpdz (5.51d) 
z + + z)(/+ + f~){t,x + z,p) dpdz 

+c-'^ J {H+{t,z)r{t,z)-H~{t,z)j+{t,z))dz + 0{c-'). 

Analogous computations lead to 
B{t,x) = 0,ptic-') (5.52a) 
\z\'''^z A pf{t, x + z,p) dp, dz + Ccpt(c"^) (5.52b) 



= c 



= c ^ 1 1 \z 

^-3 



z Ap — c ^-2 A p{z ■ pY^f{t, X + z,p) dp, dz 
z]''^ i^{z ■ p)z A {■ ■ ■ ) + {z A p) IS) z^E{t,x + z){f^ + f'){t,x + z,p) dpdz 
J ^ ^ {H^{t,x + z)j~{t,x + z) - H'{t,x + z)j^{t,x + z)),dz 
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+api(c-'). 



(5.52c) 



Now we have to prove the error estimates step by step, in fact for proving the 1.5PN 
approximation we need to know that E = Eq + Ocpt{c~^)- 

Lemma 5.1 (Newton approximation) For all < T < T and R > there are con- 
stants M{T) and M{T, R) such that 



ir(t,x,p)-/o^(t,x,p)| < M(T)c-^ (xeR'), 
\E{t,x) - Eo{t,x)\ < M{T,R)c-^ (l^^l < ^), 
\B{t,x) - c-^Bi{t,x)\ < M{T,R)c-^ {\x\ < R), 

for allpe M^ t G [0,T], and c>Pi. 

The proof will not be carried out here because with some rather obvious modifications it 
closely follows the lines of Section HI using the representations ()5.47|1 . ()5.52a|l and 



E^{t,x) = Eo{t,x) + 0{c-'^) and B^{t,x) = c~^Bi{t,x) + 0{c- 



(5.53) 



the last two estimates are an easy consequence of ()2.24j) and the bounds of Proposition 
Oand Lemma O 

Remark 5.2 // we replace (f^,Eo) with a solution of IVY)) the Newton approximation 
can be obtained in the same way. 

As a last step in the proof of Theorem 13.41 we have to provide an estimate of the 
differences of the third order terms. Using the Calderon-Zygmund inequality jU Thm. 
4.31] and ()2.22|) we have the estimate 



'H{z)pfQ^{t,- + z) dz 



cz 



L2 



|pfo)(^,-)IU2 



with a certain constant Cqz- Hence 



[H^j- - - H+j, + H-j^) (t, z) dz 

< ii(^+-/7+)(t,-)iiL-iir(t,-)iU2 + iiif-*-(t,-)iu-ii(r-jo")(t,-)iiL2 
+ll(^- - H-){t, . ||j+(t, .)u + ||if-(t, .)u . ||(j+ - jo^)(t, 

< C^Czf||(p+-Po^)(t,-)l|L-||r(t,-)l|L2 + ||p+(t,-)||L-||(r-Jo")(t,-)llL2 



+ II(P -Po)(^, 



l|j+(i,-)IU2 + ||Po(t, 



\\{J^-Jt){t. 



< CczQ\t)M{r^f{p*y 

< MQ\t) 



(5.54a) 



where 



Q°(t) :=max{|/=t(t,x,p)-/o^(t,x,p)| : a; G M^p G M^ ± G {+, -}} 
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Here we used (15 .4811 and Assumption I3.3r a) in combination with ()3.33|) and ()3.32|1 . In 
view of the formulas ()2.25a|l . (j5.51dj) . Lemma IKTT] ()5.53j) . ()5.54aj) and proceeding analo- 
gously to Section mU we conclude the estimate 

\E{t, x) - E'^it, x) I < Mr{c-^ + Q{t))) 

for all |t| < T, |x| < i? with a constant M/j independent of c where 

Q(t)=max{|/±(t,x,p)-/|(t,a;,p)| : a; G M^p G M^ ± G {+, -}} 

Employing the formulas (j2.25b|l and (j5.52c|l together with an estimate corresponding to 
fl5.54ap we conclude 

\B{t, x) - B'^it, x) I < Mr{c-' + Q{t))) 

for all |t| < T, < i? with a constant M/j independent of c. Proceeding analogously to 
Section 14.21 finishes the proof of Theorem 13.41 □ 



6 Appendix 



6.1 Representation of the approximation fields and 

Here we will derive the representation formula ()4.3fia|l for the approximate field E^ from 
fl2.24|) . Since the calculations for the electric and the magnetic field are quite similar, we 
will only analyse in detail the electric field and simply state the result for its magnetic 
counterpart. Actually, the representation of the electric field exhibits more difficulties 
than the representation of the magnetic field, due to the presence of the radiation term. 
From we recall E^ = Eq + c'^Ea + {2/3)c-^D^^\ where 



E^it^x) = — J \z\ '^z po{t,x + z) dz, (6.55a) 
E2{t,x) = -J z d^pQ(t,x + z) dz — J \z\^^dtjo{t,x + z) dz 



\z\ '^z p2{t,x + z) dz, (6.55b) 

D[3](t) = D^'^ = dtjj Eo{t,x){f+ + f^){t,x,p)dpdx (6.55c) 
cf. ( rrVPe ) and ()LVP|) . Firstly, we write 

D^^\t)=dt [[ Eo{t,x + z){f+ + U){t,x + z,p)dpdz, xeR\ 



Secondly, we split the domain of integration in {\z\ > ct} and {\z\ < ct}; note that the 
exterior part exactly gives E^^^ in ()4.36a|) . To handle the interior part {\z\ < ct}, in the 

sequel denoted by -Ej^^, we fix i? > and < T < min{T, T} and put 

= R + Mq and p^ = Mq, 
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recall that Mq = max{Mi{T), MsiT), M^iT)}; and are only depending on our basis 
constants tq, Sq, see (12 .2^1) . and of course of R. For |x| < R and < t < T, if |z| > 
or \p\ > p\ than f^(T^x + z,p) = f^{T,x + z,p) = for all < r < T by Theorem 
I2.ir b) and Lemma l2.2r a). This argument shows that we can replace _/|^|<cj/ dpdz by 
I\z\<mm{ct ri} I\p\<pi dp dz in the integrals defining the interior part. In other words, we 
may always assume that both \z\ and \p\ are bounded with a bound only depending on 
the basic constants and R, but not on c. Next we expand the densities w.r.t. t about the 
retarded time t{z) := t — c~^\z\ and obtain e.g. 



\z\<ct 



z\ '^z po{t,x + z) dz (6.56a) 
[ \z\-^z ((1 + c~Mdt + ll2c~M^d^^ + l/Qc-M^'d^t) 

J\z\<ct ^ 

Po(t>), a; + z) + |2|^0(c-^)) dz 
I \z\-^z (1 + c-^\z\dt + l/2c-'^\z\'^d1 + l/6c~39f ) pQ{i{z),x + z) dz 

J\z\<ct 



^cpt 



where Theorem 12.11( c) was utilized and hence M2(T) enters the bounds on 0{c~^) and 
Ocpt{c~'^). In the same manner we expand the terms in ()6.55b|) up to first order and 
the term from ()6.55cj) up to zeroth order, employing Theorem 12. iT c) and Lemma f2.2r b). 
therefore also M^iT) enters the bounds. 



[ {l/2zd^po-\z\-'dtjo-\z\-hp2}it,x + z)dz (6.56b) 

J\z\<ct 

/ {(^l + c-'\z\dt){l/2zd^po-\z\-'dtjo-\z\-hp2)}{i{z),x + z)dz 



2/3c"M dt{EQ{pl\ p^)\{t,x\z)dz (6.56c) 

J\z\<ct 

= 2/3C-3 / dt [Eo{p^ + p-)] (t», x + z)dz + 0,pt{c-^). 

J\z\<ct 

Next we sort the terms according to their orders in c~^. To zeroth order we have 

E^^^it, x) = - I \z\-'^zpo{t{z),x + z)dz + 0{c-^). (6.57) 

J\z\<ct 

In the first order we have the contribution —c~^J,,^^^\z\~^zdtPoii{z),x + p) dz, which 



gives, using that dtpo + V ■ jo = also holds true for (jrrVPJ) and further integration by 
parts. 



c 



^ / \z\ ^zV^ ■ jo(- • ■)dz 

' / \z\-\V, + c-'zdt) ■ [jo(t»,a; + z)] dz 

J\z\<ct 



'\z\ 
C 
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c ^(ct) ^ / zz ■ jo{0,x + z)ds{z) - c ^ / \z\ ^{-2zz ■ jo + jo){. . .) dz 

\z\=ct J\z\<ct 



+C-M \z\^^zz-dtjoi.--)dz, (6.58a) 

J|z|<ct 

where (...) = (i(z),x + z). Hence, the first order term of can be written as 

c-\ct)-^ zz-jo{0,x + z)ds{z)-c-^ \z\'^{-2zz-jo+3o)i---)dz (6.58b) 

J\z\=ct J\z\<ct 

To continue we collect the terms of second order from ()6.56b|l . ()6.56a|l . ()6.58aj) . 



/ (1/22:9^^0 - \z\ ^dtJo - \z\ ^zp2 - l/2zd'^po + \z\ hz ■ dtjo) (• • •) dz 

JlzKct ^ ^ 



'|z|<ct 

Since the terms containing second derivatives cancel each other we start with the first 



order time derivatives. Utilizing (jrrVPJ) and integration by parts we calculate 



c / \z\ {zz ■ dtjo - dtjo){- ■ •) dz 

'\z\<ct 



[ \z\-\z ®z- l){Eo + c-32/3D[3])^^+ ^ p--^^^ ^ ) 

J\z\<ct 

+c~^ / / \z\^'^fo{...,p)[ — 3z{z-py + 2z-pp + zp'^)dpdz 

J\z\<ct J 

-c-(.r'/ /(..p)(...p-p)/„"(. + ..p)*d. 



z\=ct 



—c^ / \z\ ^{z ■ p){zz ■ p — p)dtfo{- ■ ■ -iP) dpdz. (6.59a) 

J\z\<ct J 

Note that since the bounds in Theorem l2.ir b)(c) also imply that 2/3c~^ J\z\<ct \ 

1)D[^](Pq + Po )(• ■ ■) dz = Ocpt{c~^) this term can be dropped. Using the same argument 

all terms containing appearing in the sequel are at least Ocpt{c~^) and hence will be 

dropped without comment. Thus, the contribution of second order in is 
/ \z\-\z®-z-l)Eo{...)U^ + fo){....p)dpdz 

J\z\<ct 

J\z\<ct L ^ 

-c-Hctr f ■ rt(.-^= ■ P - P)fS(^ + ..P) ipi^. (6.59b) 

J\z\=ct J 

We continue by collecting the terms of third order from (j6.56bj) . ()6.56cj) . ()6.56a|l and 
[ {(-1/6 + l/2)zdfpo - d^jo - \z\-'zdtP2 + 2/3dt[Eo{pt + Po ] )(• • •) dz 

J\z\<ct ^ 

z\~^{z ■ p) {zz ■ p)dtfo{- ■ •) dpdz (6.60a) 



\z\<ct 
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Using JrrVPel ) and integration by parts we compute 

l/3c-^ / 2<9fpo(- • •) dz (6.60b) 

J\z\<ct 

= l/3c-^ / dfjo{. ..)dz- l/3te~2 / zz ■ d'fjo{0, x + z) ds{z) + 0(0"^). 

J\z\<ct J\z\=ct 

Summing the terms with two time derivatives from ()6.60a|) and ()6.60b|) we get 
-2/3C-3 / d^jo{...)dz = -2/3c-' [ dt[{Eo + {2/3)c-'D^'^){pt + p,)]{...)dz 

J\z\<ct J\z\<ct 

(6.60c) 

+2/3c-^ / I z ■ ppdtfo{0,x + z,p)dpds{z) + 0{c-^). 

J\z\=ct J 

Now we collect the first time derivatives from ()6.60ap . Starting with the contributions 
containing /q we have 



z\<ct 
-3 



z\<ct 



\z\ ^{z{z-py - {z ■ p)p)dtfo{- ■ ■ ,p)dpdz 
\z\^'^ (2{z ■ p)z z — z ■ p — p 1^ z) 

[Eo{. . .) + (2/3)c-3D[3](t»)) (/o+ + /o-)(. ..,p)dpdz 
—c^'^ / / \z\^^ fo{. . . ,p)(^ — Az{z ■ p)^ + 3p{z ■ p)"^ + 2z{z ■ p)p'^ — pp"^) dpdz 

J\z\<ct J 

+c~^{ct)-^ 1 I {z{z -pf -p{z ■pf)f°{x + z,p)dpds{z) + 0{c~^) (6.60d) 

J \z\ =ct J 

and secondly, for the terms containing /2 we compute employing ()LVP|) 
-c"^ / \z\~'^zdtp2{- ■ •) dz 

J\z\<ct 

= -c'^ I \z\-'^f2{---,p){-2z{z-p)+p)dpdz 

J\z\<ct J 



+c-^ctr I I 

J\z\=ct J 



z{z ■ p)f2{x + z,p) dpds{z) 



+1/2C-3 [ [ \z\-^fo{...,p){-2z{z-p)p^ + p^p)dpdz 

J\z\<ct J 

-l/2c'^{ct)-^ I z{z-p)p^f^{x + z,p)dpds{z) + 0{c-^). (6.60e) 

J\z\=ct J 

Summarizing ()6.60a|) - (j6.60e|) we can identify the contribution in the third order of E^^. 

/ / { - l/3tc-'^z{z ■ p)d1 + 2/3c-^{z ■ p)pdt 

J\z\=ct J 

+c~^{ct)~'^(^z{z ■ pf -p{z-pY - l/2z{z ■ p)p^^ }/o(0,x + z,p) dpds{z) 
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.-3 



+c"^(ct) ^ / / z{z ■ p)f2{x + z,p) dpds{z) 

J\z\<ctJ ^ ^ 

W / \z\-'{2{z-p)z^z-z-p-p® z}Eo{...)Uo + fo ){■■■, P)dpdz 

J\z\<ct J 

+c'^ I \z\-^f2{...,p){2z{z-p)-p}dpdz. (6.60f) 

JlzKct J 



Therefore, if we use fl6.57p . ()6.58b|) . ()6.59b|) . ()6.60f|l and add some terms of order 0{c ^ 
containing as e.g. ij^|<ct / kl"^^^^ 
composed as it is claimed in ()4.H(ja|l 
Similar calculations for yield 



z\ ^z(z ■ pY f2{- ■ ■) dpdz it turns out that can be de- 



(6.61) 



where 



z\>ct 



z\~^z X (jo + - l/2c~^zd^jo) {t, x + z)dz 



z\<ct 



\z\ ^Li{z,c ^p)fR(i{z),x + z,p)dpdz 



/ / \z\ 'L2iz,c-^p)Eoif^ + f^){i{z),x + z,p)dpdz 

J\z\<ct J 



(ct) / / Ls{z,c p)fQ{x + z,p)dpds{z) 

'\z\=ct ■ 



and the kernels are given by 

Li{z,p) = z X p — 2z X p(z ■ p)3/2z X pp^ + 3z X p(z ■ p)'^ 

Lz{z,p) 



z X [■■■)— z ■ pz X {■■■) — z ■ pz ■{■■■ ) G 
—z X pz ■ p + z X p{z ■ py. 



6.2 Representation of the Maxwell fields E and B 



In this section we will verify the representation formula ()4.38a|) for the full Maxwell field 
E by expanding the respective expressions from |13t 1^ to higher orders. Once again the 
computation for the corresponding magnetic field B is very similar and therefore omitted. 
Let {f,E,B) be a C^-solution of JRVMeD with initial data according to 

dEl- We recall the following representation from ^ (A13), (A14), (A3)], 



E = Ed + Edt + Et + Es, 



(6.62) 



where 



Eoit^x) = dti — E°(x + ctuj)duj] +— dtE(0, x + ctuj) du, 
V4vr J\^\=i J 47r 

EDT{t,x) = {cty^ / / KDT{z,p)fo{x + z,p)dpds{z), 

J\z\=ct J 
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Erit^x) = / KT{z,p)f{...,p)dpdz, 

J\z\<ct J 

Es{t,x) = / \z\-' [ Ks{z,p){E + c-'px B)i...){f + + /-){..., p)dpdz 

J\z\<ct J 

The kernels are given by 

Kdt{z,p) = -{1 + c'^z ■ py^{z - c''^{z ■ p)p), 
Kt{z,p) = -{l + c-h-p)-^{l-c~^p'^){z + c-^p), (6.63a) 

Ks{z,p) = -{l + c~^z-py^{l + c~'^p'^)-^^^ {l + c'^z-p) + c~^{{z-p)z~p)®p 

-{z + c-^p)^z eR^""^ (6.63b) 

Next we expand these kernels in powers of c~^. According to ()3.27ap the p-support of 
/^(t, X, ■) is uniformly bounded in x G and t G [0,T] by M5. Thus, we may suppose 
that \p\ < M^{T) in each of the p-integrals. Therefore as long as /^(. . . ,p) 7^ we have 

Kt{z,p) = K,{z,c-'p) + 0{c-^) 
Ks{z,p) = K2{z,c-^p) + 0{c-^) 
Kdt{z,p) = -z + K:,{z,c~^p) + 0{c~^). 

Recall that /°'^ = if |x| > ro, see fl2.23j) and utilizing this as well as 



-{Ct)-^ / / 0{c-^)lB(Q,r,){x + Z) dpds{z) 

J\z\=ct V|p|<pt 

= (ct I lB(,0,r,){x + Ctuj)ds{uj)\O{c-^) = O{c-'') 
V J\u>\=l J 

by Lemma 1], uniformly in x G M^, t G [0, T], and c > 1, we arrive at 
EoTit^x) = (ct)-' [ [ {-z + K,{z,c-'p)){f° + c-^f°){x + z,p)dpds{z) + 0{ 

J\z\=ct J 

(6.64) 

Concerning Et, we note that f^{t,x,p) = for |x| > M^(T) and t < T, see fl3.27ap . 
Since, by distinguishing the — y| > 1 and \x — y\ < 1, 



c 



\z\ 1b{o,M5{t)){x + z)dz = \x-y\ lB{o,M5{T)){y) dy = 

z\<ct J \x—y\<ct 

uniformly in x G M^, t G [0,T], and c > 1, similar computations as before show that 

Erit^x) = / \z\-^ Ki{z,c^^p)f{...,p)dpdz + 0{c^^). (6.65a) 

J\z\<ct J 

In the same manner, elementary calculations using also ()3.27b|) can be carried out to get 
Es{t,x) = c-^ [ \z\'' [ K2{z,c-'p)E{. . .){f+ + f-){. . . ,p) dpdz + 0{c-^). (6.65b) 

J\z\<ct J 

Observe that we have proved 

ET + Es = E,r,t + Oic-^), (6.65c) 
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see fl4.38d|l . Next we consider the data term 

Eoit, x) = dt(— / E°{x + ctu) duj) + — / 9(^(0, x + ctu) du = I + 11 (6.66) 
47r i|^|=i 47r J\^\^^ 

By Maxwell's equations we have dtE{^,x) = cV x B°{x) — 47rj(0,x). Recall that B° = 
c-^Bi + c-^B^ + c-^Bcjree- Using A = -V X V X +VV- and dtpo + V ■ jo = as well 
as dtP2 + V ■ j2 = 0, the latter an easy consequence of ()LVP|1 . it is easy to check that 
V X i?i = OiEq + 47rjo and V x i?3 = dtE2 + 47rj2- Employing the definition of Eq, E2 
and ()6.7()aj) - ()6.7()c|l below we compute 

— / dtEo{0,x + ctuj) duj = —— / / \z\^'^zdtpo{0, X + ctuj + z) dz duj 

4^'^ J\cu\ = l J\u,\ = lJ 



t 



~ J dtpoiO,y) J \y ~ X - ctuj\'^{y - X - ctuj)dujdy 
-t I \z\^'^zdtpo{^.,x + z) dz (6.67a) 

J\z\>ct 

dtE2{0,x + ctuj)duj (6.67b) 

Ul\=l 

^ I I {l/2zdlpQ-\z\~^d'liQ-\z\~'^zdtp2]{'^,x + ctuj + z)dzduo 
^ / zdfpo{0,x + z)dz+^ I zdfpo{0,x + z)dz 

C J\z\<ct -^C y|2|>ct 

I \z\^'^zdfpo{0, x + z)dz - ^ I d^jo{0, x + z) dz 

J\z\>ct C J\z\<ct 

/ \z\^^d'^jo{0,x + z) dz - / \z\^^zdtP2{0,x + z) dz. 

C J\z\>ct C J\z\>ct 



Using JrrVPcD we also have 

/ zdlpo{Q, x + z)dz-\ I d^joiO, x + z)dz (6.67c) 

J\z\<ct C J\z\<ct 



~ I zz ■ d^joiO, x + z)dz- ^ I dtjo{0, x + z) dz 

/ ZZ ■ dtjo{0,x + z)dz + ^ I nz-p)dtjo{0,x + z)dz 

Ac J\z\=ct Ac J\z\=ctJ 

2 



2 <ct 



dt [(Eo + (2/3)c-3dP1) (Po + Po )] (0, x + z) dz. 



Furthermore, according to ()2.23p and ()3.26p and the definition of jo and j2 the support 
of j := (j — jo — c~^j2)(0, •) is bounded by Tq and 

m = [ \{p-p+ |^)/o° + c-'{p - p)f^ + c-'pf^j] {x,p) dp = 0{c-'). 

J \p\<ro L J 

Utilizing |24[ Lemma 1] we conclude 

^ / A7ij{0,x + ctiu)duj = O{c-^). (6.67d) 
471- J\w\=l 

32 



Now we return to the contributions coming from I in ()6.66j) . Note that in p304/305,(5.21)] 
the contributions coming from Eq and E2 are already determined to 

\z\^'^z po{0,x + z) dz + - c"'^ I z d^po{0,x + z) dz 

z\>ct ^ J\z\>ct 



— I \z\ "^z Pq{0,x + z) dz — c ^ j \z\ ^dtjo{0,x + z) dz 

2 J\z\>ct J\z\>ct 

+ {ct)-^ z(po + c-V2)(0,x + 2)ds(z). (6.68a) 

J\z\=ct 

Concerning the two remaining terms we have 

^sl^^'^^^^^sl / dt[E,{pt + pM^,z)dz (6.68b) 
= / 9t[Eo(p++Po)](0,x + 2;)dz + -^ /" dt[Eo{pt ^ p^)\{Q,x + z)dz, 

<JC J\z\>ct 'JC J\z\<ct 

note that the last term in ()6.68bj) cancels the last term in ()6.67c|l : and 



d,{^ f c-^'E^j^Jx + ctu) duj) (6.68c) 



^-4 r ct f 

4vr J\^\=i '■' J\^\=i 



according to ()3.26|) . Lemma 1]. Thus, combining ()6.64|) and ()6.66p - ()6.68c|) the repre- 
sentation formulas (j4.38bj) and (j4.38djl are proved. 

6.2.1 Representation of B 

We simply state the representation of B. 

B = Be.t + Bint + Bbound + 0{c-^) (6.69) 

Be.t{t, x) = I {\z\-H^ {f - tdtjo - tV2jo) 

J\z\<ct 

-l/2rt X d'fjo}{0,x + z) dz 
Bint{t,x) = / \z\~^Li{z,c'^p)f{i{z),x + z,p)dpdz 

J\z\<ct J 

[ [ \z\-^L2{z,c-^p)E{f^ + f-){i{z),x + z,p)dpdz 

J\z\<ct J 
Bbound{t-i x) = -Bfe^und 

6.3 Some explicit integrals and a lemma 

We point out some formulas that have been used in the previous sections. For 2; G 
and r > an elementary calculation yields 

I 1-1 r f 4:7^^^ : r > ,„ „ . 

\z — ruj\ duj = < , : : , ^ • (6.7(Ja) 

' ' I Anlzt^ : r < \z\ ^ ^ 
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Differentiation w.r.t. z gives 



OJ =1 



\z — ruj\ "^(z — tlS) duj = \ , I I 9_ ' . (6.70b) 



Similarly, 

/ \z ~ ruj\du) 

J\i0\ = l 

and thus by differentiation 



Anr + ^z'^r ^ : r > |2;| 



Att\z\ + ^r^\z\~^ : r < \z' ' 



^"^ z : r > l^;! 



z — ruj\ ^ {z — ruj) duj = < _ ^ . (6.70c) 



'\ui\=i I 4:7rz — ^r^\z\ ^z : r < 

Finally, for ^ G M''^ \ {0} also 

j \z -v\-^\v\-^vdv = 2txz (6.70d) 

can be computed. 



Acknowledgments: The author is indebted to M. Kunze, G. Rein, A. Rendall and 
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